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Abstract 

We demonstrate the effect of single-electron tunneling (SET) through a carbon nanotube quan- 
tum dot on its nanomechanical motion. We find that the frequency response and the dissipation of 
the nanoelectromechanical system (NEMS) to SET strongly depends on the electronic environment 
of the quantum dot, in particular on the total dot capacitance and the tunnel coupling to the metal 
contacts. Our findings suggest that one could achieve quality factors of 10 6 or higher by choosing 
appropriate gate dielectrics and/or by improving the tunnel coupling to the leads. 



Carbon nanotubes (CNT) have become an essential building block for nanoelectrome- 
chanical systems (NEMS). Their low mass and high Young's modulus allow for instance 
ultrasensitive mass [TH5] or force detection jU ]5\ (electric and magnetic) over a wide range 
of frequencies and its small diameter enables even single- molecule detection [SHI]. More- 
over, CNT devices exhibit remarkable electronic transport properties, ranging from Kondo 
physics [8 J to Coulomb blockade at high temperature (9] . It was recently demonstrated that 
a CNT NEMS' nanomechanical motion at very low temperature (i.e. in Coulomb block- 
ade regime) is strongly affected by the electronic transport through the CNT quantum dot 
(QD), and vice versa: For instance, single-electron tunneling (SET) caused a frequency soft- 
ening and increased dissipation when tuning the CNT dot's potential through a Coulomb 
peak PEHE]. 

Here we demonstrate that the response and dissipation of a CNT NEMS at very low 
temperature induced by zero bias SET through the CNT NEMS-QD critically depends on 
the dot capacitance, the tunnel coupling to the metal leads and temperature. We studied 
the frequency and dissipation response of nanomechanical bending modes to zero bias SET 
in suspended CNT devices with tunable tunnel couplings and different gate dielectrics, i.e. 
different dot capacitance. We observe that SET causes a frequency softening for small 
dot capacitance and/or tunnel coupling, whereas a frequency hardening or no frequency 
modulation is observed for large dot capacitance and/or tunnel coupling. We show that 
the dissipation of the CNT NEMS is mainly dominated by the capacitance, when electron 
tunneling through the dot is suppressed (i.e. in the Coulomb valley), whereas it is limited 
by the mean tunneling rate T and the conductance, when electron tunneling through the 
dot is allowed (i.e. on a Coulomb peak). Finally we demonstrate that the tunnel current is 
the dominant dissipation mechanism in CNT NEMS at low temperature. Our findings are 
in fair agreement with a theoretical model provided previously [101 EH] ■ 

Low capacitance CNT quantum dots (Cdot ~ 20 - 40 aF) were obtained by using silicon 
dioxide as gate dielectric [Fig. 1(a)]. First, source-drain electrodes are patterned by optical 
DUV lithography and e-beam evaporation of Mo (20 nm) and Pt (160 nm) on 500 nm 
of thermal Si02- To ensure the suspension of the CNT, 150 nm of Si02 are dry etched in 
CHF3 plasma. High capacitance CNT quantum dots (Cdot ~ 160 - 260 aF) are obtained with 
high-K gate dielectrics Al 2 03 [Fig. 1(b)]. First, a 1 yum-wide metallic local gate is patterned 
by optical DUV lithography and e-beam evaporation of Mo (20 nm) on 300 nm of thermal 
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FIG. 1. False color SEM image of (a) a low capacitance device based on a SiC>2 covered Si 
backgate (grey color) and (b) a high capacitance device based on an AI2O3 covered local metallic 
gate (blue). The RF actuation signal is injected into the local metallic gate and Si backgate, 
respectively, through a home-built bias-T. As the induced mechanical motion changes the charge 
flow through the CNT quantum dot and vice versa, we can detect the CNT resonance through 
a change in zero bias conductance, (c) Mechanical resonance of a typical CNT NEMS at low 
driving power Prf = — 100 dBm. The resonance width 5f = 2.5 kHz leads to a quality factor of 
Q 140000. 



Si02- A layer of 100 nm of AI2O3 is then deposited by atomic layer deposition. Using optical 
DUV lithography and e-beam evaporation of Mo (20 nm) and Pt (160 nm), source-drain 
electrodes are aligned above the local gate. Suspended CNT's are finally grown by chemical 
vapor deposition at 800°C from a CH 4 feedstock and Fe/Mo catalyst spots patterned on the 
source-drain electrodes next to the junction. The CNT device length is approximately 800 
nm (± 50 nm) and the dot capacitances are deduced from the dot's charging energy at low 
temperature. The spread in capacitance values of nominally identically fabricated devices 
is due to variations in the nanotubes length and slack (see Supporting Information). 

The measurements are carried out in a 3 He/ 4 He dilution refrigerator with a base tem- 
perature of 30 mK. The NEMS actuation and detection scheme used in our experiment is 
similar to the one used by Steele and co-workers [TT1 [12]. The RF actuation signal is in- 
jected into the gate electrode via a home-built bias T. As the induced mechanical motion 
changes the charge flow through the CNT quantum dot and vice versa, we can detect the 
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CNT resonance through a change in zero bias conductance. The actuation power is kept to 
a minimum (Prf ~ -100 dBm) in order to ensure a quasi-linear regime of the oscillator and 
a high Q lorentzian resonance shape [Fig. 1(c)]. All measurements were done under zero 
bias with a standard lock-in technique. 

SET in a CNT NEMS-QD can be considered as an external perturbation to the CNT 
mechanical motion, and vice versa. This perturbation can be described as an electrodynamic 
force acting on the CNT. The contribution of this force, which is in phase with the mechanical 
motion, is responsible for a frequency modulation Af, whereas a contribution of the force 
being out of phase with mechanical motion induces a modification of the dissipation and the 
quality factor Q. Af and Q can be expressed as follows [10] 



Where T is the mean tunnel coupling of the dot to the leads, yielding T = (T d + T s ) /2, 
and Cdot is the total dot capacitance. We can estimate the tunnel coupling T for each device 
from the line shape of the Coulomb blockade peaks (see Supporting Information). Eqs. (1) 
and (2) are valid in the regime of Coulomb blockade and if T ^ f . 



dot capacitances C^ot- We studied the frequency response for comparable Coulomb peaks 
height, i.e. for comparable conductance G (V g ) ~ 20 /j,S. The tunnel coupling are on the 
order of 10 GHz. As depicted in Fig. 2(a), we observe a strong frequency softening for small 
dot capacitances, whereas for high dot capacitances a frequency hardening or no modulation 
is observed. As shown in Fig. 2(b), the experimental data are in fair agreement with the 
model. A fit with Eq. 1 yields fitting parameters C'^/k = 8 • 10" 22 F 2 /Nm and r 10 GHz, 
which are in rather good agreement with previous experiments [TU1 HU [T5] . Fig. 2(c) shows 
the (scaled) frequency modulation for different Coulomb peaks and tunneling rates T on a 
given device (device 1). It was shown previously that one can tune a CNT quantum dot 
from a Coulomb blockade regime to a regime of strong tunnel coupling by simply changing 
the CNT's gate potential [HJ [15]. Hence, we can change the tunnel coupling T by tuning 
the gate voltage of our CNT device and estimate the values of T from the Coulomb peak line 




(2) 



(1) 



First, we compare the frequency response Af (scaled with /• V g 2 ) of devices with different 
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FIG. 2. Modulation of the resonance frequency due to SET: (a) Frequency shift A/ (scaled with 
/ • Vg) for different dot capacitance Cd t> with comparable tunneling rates T. inset shows a typical 
Coulomb peak centered at AV^ = 0; (b) Maximum scaled frequency shift at maximum tunnel 
current G max (AV g = 0) vs. the dot capacitance for all measured devices. The solid red line 
corresponds to the fit using eq. 1 with C^/k = 8 • 10" 22 F 2 /Nm and T d ps 10 GHz. The data 
are in fair agreement with the model, (c) Scaled frequency shift for different Coulomb peaks i.e. 
tunneling rates T, with a fixed dot capacitance Cdot = 260 aF (device 1). The tunneling rate T is 
estimated from the lineshape of the respective Coulomb peaks (see Supporting Information) . (d) 
Maximum frequency shift A/ max at maximum tunnel current G max (AV^ = 0) vs. temperature for 
a device yielding Cdot = 160 aF (device 2). The inset shows the evolution of G max as a function of 
temperature. 



shape (see Supporting Information). The dot capacitance Cdot = 260 aF remain unchanged 
by tuning the gate voltage. As T increases, the CNT becomes stiffer, resulting in a sign 
change of Af and the transition from a frequency softening to frequency hardening, in 
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FIG. 3. Modulation of the Q-factor (scaled with / • V g 2 ) due to SET: (a) For different dot capaci- 
tances, with comparable tunneling rates V . The solid lines are guide to the eye and the inset shows 
a representative Coulomb peak centered at Alg = 0. (b) Maximum quality factor Q max at minimal 
tunnel current G (| AV g \ 3> 0) vs. the dot's capacitance for all measured devices. The solid red line 
represents the fit with eq. 2 with C g 2 /k = 1 • 10" 21 F 2 /Nm and T ~ 10 GHz. Our findings are in 
good agreement with the model, (c) Modulation of the Q-factor (scaled with / • V g 2 ) for different 
tunneling rates T, with a fixed dot capacitance Cdot = 160 aF (device 2). The solid lines are guide 
to the eye and the tunnel coupling V is estimated from the lineshape of the respective Coulomb 
peaks (see Supporting Information), (d) Maximum quality factor Q max at minimal tunnel current 
GdAV^I 3> 0) vs. temperature for a device yielding Cdot = 160 aF (device 2). The inset shows 
the evolution of G m i n function of temperature. 



qualitative agreement with Eq. 1 [Fig. 2(c)]. 

Fig. 3(a) depicts the effect of SET on the dissipation of the CNT NEMS, i.e. the Q 
factor (scaled with / • V g 2 ), for devices with different dot capacitances. From Eq. 2, we 
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expect an increase of the Q factor with increasing capacitance in the limit of suppressed 
SET through the dot. Indeed, in a region of suppressed SET (|AVg| 0),the Q factor is 
larger for devices with higher dot capacitance [Fig. 3(a)]. The tunnel coupling r ~ 10 GHz 
and the conductance G (V g ) ~ 18 /xS are comparable for both traces. Fig. 3(b) shows the Q 
factor (scaled with / • Vg) in the limit of suppressed SET (|AV^| ^> 0) as a function of the 
dot capacitance for all measured devices. The fit with Eq. 2 yields C^/k = 1 • 10" 21 F 2 /Nm 
and r « 10 GHz which is consistent with the fitting parameters of Fig. 2(b) and previous 
experiments [TUJ [III H3]. Despite the good agreement between theory and experiment, we 
observe deviations from the model for some devices in Fig. 2(b) and 3(b), which we attribute 
to variations in C g 2 jk (see Supporting Information). Moreover, we expect from Eq. 2 and 
previous experiments [TQj, [TT] an enhanced electromechanical dissipation in a given device 
if the conductance G (i.e. the tunnel current) and the tunnel resistance at the nanotube- 
electrode interface R ~ 1/T d are increased, in analogy to the current dissipation in a simple 
resistance. Indeed we observe increased dissipation (decreased Q factor) when tuning the 
gate voltage through a Coulomb peak [Fig. 3(a) and 3(c)]. In Fig 3(c) (device 2), this effect 
becomes more pronouced as we move to gate voltage regions with a smaller tunneling rate 
T, i.e. larger tunnel resistance at the nanotube-electrode interface. The dot capacitance 
Cdot yields 160 aF for this device and the height of the Coulomb peak are comparable for all 
traces (G ~ 7 fiS). Therefore, the CNT capacitance is the limiting factor for the dissipation 
in regions of suppressed SET, whereas the mean tunnel coupling and the tunneling current 
itself define the dissipation in regions of strong SET. We conclude that the electron transport 
through the CNT is the main dissipation mechanism in our nanoelectromechanical system. 

Finally we study the temperature dependance of the frequency response and the dissipa- 
tion. Fig. 2(d) depicts the frequency shift A/ max for G max = G (AVg = 0) whereas Fig. 3(d) 
shows the quality factor Qmax for G m i n = G (| AVg| ^> 0) as function of temperature. A fit of 
the data yields A/ max ~ T and Q max ~ \/T. As temperature increases the Coulomb block- 
ade peak broadens and becomes smaller (see supporting information Fig. S3). As a result 
C max = G (AVg = 0) decreases whereas G min = G (|AVg| ^> 0) increases with temperature, 
as depicted in the insets of Figs. 2(d) and 3(d). From the data fit we obtain G max ~ — T and 
Gmi Q ~ T. It was previously demonstrated that carbon nanotubes show Luttinger liquid 
behavior [TBTfPB"] . where the conductance follows a characteristic power law dependence as 
a function of temperature G(T) ~ T a , a being related to the Luttinger parameter g by 
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a = (g — 1 + g — 2)/2 [TBI ELY]- We ^ n( ^ a ~ 1 anc ^ a Luttinger parameter of g = 0, 1, which 
is close to previous measurements |16] and theoretical predictions [TT[ HE] • Comparing the 
temperature dependance of the conductance, the quality factor and the frequency shift we 
finally obtain A/ max ~ -G max and Q 

max ^ 1 f Gyr\\ j\ j which is in full agreement with Eqs. 1 
and 2. We can conclude that the current is the dominant dissipation mechanism for carbon 
nanotube based NEMS in a Coulomb blockade regime at cryogenic temperatures. 

We have demonstrated that the response and the dissipation to single-electron tunneling 
of carbon nanotube NEMS at low temperature depends on the dot capacitance and the 
tunnel coupling to the leads, the tunnel current itself being the main dissipation mechanism 
in the system. By choosing high-K dielectrics (Hf02 or Zr02) and improving the tunnel 
contact to the metal leads, one could significantly enhance the quality factor Q to values 
exceeding 10 6 . It was proposed in theoretical calculation that one can use such high— Q 
CNT NEMS as magnetic torque or force detectors for nanoparticles [5J or single molecule 
magnets [T9l [20] grafted to the CNT NEMS. In principle one can achieve a sensitivity of one 
(Ib at low temperature, whereas the best magnetometers, for instance the micro-SQUID, 
only have a sensitivity of 10 3 /xb [21] • Finally, such molecular quantum spintronic device 
would allow single spin manipulation on a molecular level. 
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1. Estimating quantum dot capacitances of CNT NEMS 

The total dot capacitance is obtained from bias spectroscopy measurements. The total 
dot capacitance Cdot is related to the charging energy E c of the dot via Cdot = dE c . The 
charging energy is deduced from the size of the Coulomb diamonds (fig. SI a). Alternatively 
one can determine the capacitances C s and Ca from the slope of the diamond edges and the 
gate capacitance C g from the distance between 2 degeneracy points (Fig Sib). The total dot 
capacitance then yields Cdot = C s + Cd+ C g . 

For SiC>2 based devices we obtain capacitance between 20 and 40 aF, whereas for 
AI2O3 based devices the values range from 160 aF to 260 aF. The spread in capacitance values 
for nominally identically fabricated devices can mainly be attributed to different device 
lengths (ranging from 700 nm to 850 nm) and therefore to variations in charging energies. 
Electron microscopy revealed slack on most devices, suggesting a low built-in tension which 
is characteristic for CVD grown CNT 's. 1,2 As the slack may be a different for each device, the 
distance between the CNT and the gate can vary from one device to another, which would 
also contribute to the spread in capacitance values. 
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Fig. SI: Estimating the dot capacitance via bias spectroscopy: (a)(c) Coulomb diamonds 
for two AI2O3 based devices with a length of 780 nm and 820 nm, respectively. (b)(d) 
Coulomb diamonds for two SiC"2 based devices. 



2. Estimating the tunnel coupling T in CNT quantum dot 

In the regime of Coulomb blockade, the energy level broadening in the carbon 
nanotube quantum dot is given by the electronic temperature and the tunnel coupling to the 
leads. The electronic temperature in our system is around 150 mK, therefore at least an order 
of magnitude smaller than the expected tunnel couplings in carbon nanotube junctions I>> 
kT. In this limit, the Coulomb peak line shape can approximated as follows 3 : 

^ 2 ) 2 (SI) 



G max a 2 (V g -V gfi ) +(r/2) 2 

Where G max is the maximum of the Coulomb peak, T = (T s + Td)/2 the mean tunneling rate 
through the quantum dot, F g ,othe center of the peak and a an energy scaling factor given by: 

a = = 1 (S2) 

1.36 dV g \36C S +C d +C g 

We used Eq. SI to fit the Coulomb peaks and estimate the mean tunnel coupling T for devices 
1 and 2 studied in Fig. 2 (c) and Fig. 3 (c) of the main text. The results are depicted in Fig. S2 
with the same color code than in the main text figures. For device 1, studied with respect to 
the frequency shift Af, we can tune the tunnel coupling T from 30 p,eV to 55 p,eV (7.25 GHz 
to 13.3 GHz) [Fig. S2 (a)], leading to a stiffening of the CNT [Fig. 2(c)]. For device 2, studied 
with respect the Q factor modulation, we can change T from 20 p,eV to 90 |j,eV (4.8 GHz to 
21.5 GHz), leading to an increase of the quality factor [Fig. 3(c)]. 

Although Eq. (1) is only an approximation (with an error estimated to 7%), it is 
sufficiently accurate to show a dependence of the frequency shift and the dissipation on the 
tunnel coupling in qualitative agreement with Eq. 1 and 2 of the main text. 
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AV (mV) 

Fig. S2: Estimating the tunnel coupling T in the limit of tunnel broadened Coulomb 
blockade, a) Device 1 studied in Fig. 2(c) of the main text with respect to the frequency shift 
Af. b) Device 2 studied in Fig. 3(c) of the main text with respect to the quality factor Q. The 
color code is identical to the one of the main text figures. The estimation is accurate enough to 
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demonstrate a qualitative agreement of our data in Fig. 2 (c) & 3 (c) with the model and Eq. 1 
& 2 of the main text. 



3. Accuracy of fitting parameter C g Ik 

In Fig. 2b and 3b of main text, we fitted the experimental data with equation 1 and 2, using r 
and C g ' 2 /k as fitting parameters. We obtain average value for C g ' 2 /k = 10" 21 F 2 /Nm and P= 10 
GHz, which is in fair agreement with previous experiments. 4 ' 5 Although the model fits the 
experimental data with good accuracy, we can observe certain deviation from the model for 
some devices, which could indeed be attributed to variations in C g /k: 

2 2 2 

According to ref. 4, C g /k is given by C g ' /k = (C g /Z) /k in good approximation, with C g being 
the gate capacitance, Z the CNT-gate distance and k the CNT's spring constant. 
We studied two different device layouts: SiCh based CNT QD with a gate capacitance C g = 
2±1 aF and Z = 700 nm and A1 2 3 based CNT QD with C g = 3±1 aF and Z =300 nm. C g is 
estimated from the voltage separation of two adjacent Coulomb peaks (see Fig. SI) and Z is 
given by the sum of the oxide thickness and thickness of the metal electrodes. With an 
average spring constant k = 10" 3 N/m (value extracted from ref. 4 and 5) we obtain values for 
C g ' 2 /k = (6±2)T0~ 21 F 2 /Nm, which is close to values obtained from the fit in Fig. 2b and 3b. 

The above estimation shows that we indeed have variations in C g ' 2 /k, caused by a small spread 
in capacitance values C g and CNT-gate distance Z, which could ultimately lead to deviations 
of some devices from the model fit. 
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4. Power law temperature dependence of the zero bias conductance 

As temperature increases the Coulomb blockade peak broadens and becomes smaller 
(see Fig. S2). As a result G max = G(AV g = 0) decreases whereas G m ; n = G(\AV g \ = -1 mV) 
increases with temperature, which is depicted in the insets of Figs. 2(d) and 3(d) of the main 
text. Both follow a power law dependence, for instance G max ~ -T and G m ; n ~ T, as expected 
for Luttinger liquids like carbon nanotubes. 6 ' 7 ' 8 
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Fig. S3: Evolution of the zero-bias conductance as a function of temperature (from 20 
mK to 1 K) for the device studied in Fig. 2(d) and 3(d) of the main text. 
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